THE BELL SYSTEM 
TECHNICAL JOURNAL 



DEVOTED TO THE SCIENTIFIC AND ENGINEERING 
ASPECTS OF ELECTRICAL COMMUNICATION 



Volume 60 May-June 1981 Number 5 

Copyright © 1981 American Telephone and Telegraph Company. Printed in U.S.A. 

Dispersionless Single-Mode Lightguides 
With a Index Profiles 

By U. C. PAEK, G. E. PETERSON, and A. CARNEVALE 
(Manuscript received November 1 4, 1 980) 

By means of a numerical solution to Maxwell's equations, we 
calculate those parameters necessary to design and fabricate single- 
mode lightguides. These include optimum core radius and profile 
parameter a. In the design the dispersion is minimized by varying 
the core radius while the relative index difference A, the wavelength 
X, and profile parameter are held fixed. A comparison of calculated 
dispersion with experimental data shows excellent agreement. 

I. INTRODUCTION 

Recently, single-mode lightguides have been developed that achieve 
transmission losses as low as 0.5 dB/km and 0.2 dB/km at wavelengths 
of 1.3 /mi and 1.55 jum, respectively. 1 When the total lightguide disper- 
sion is reduced to zero at the operating wavelength, a transmission 
system can be realized with wide repeater spacings and extremely large 
bandwidths. 2 " 4 In fact, band widths in excess of 1 GHz/ 100 km are 
expected. Clearly, such a lightguide system would be ideal for undersea 
cable and other long-distance transmission applications. 

The design of single-mode lightguides with zero total dispersion 
requires an accurate description of this parameter in terms of profile 
shape, materials properties, and core radius. This necessitates a solu- 
tion of Maxwell's equations for the fundamental HEu lightguide mode 
as well as the TE \ and TM 0i modes. 

For radially inhomogeneous media, it is usually not possible to 
obtain these solutions as analytical expressions of a closed form. To a 
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very large extent the design of single-mode lightguides has in the past 
focused on rectangular-shaped index profiles. This is because Max- 
well's equations are very much easier to solve in this case. When 
nonrectangular profiles have been considered, numerous approxima- 
tions have often been employed. 

To avoid these difficulties, we have developed a numerical technique 
to obtain exact solutions to the vector form of Maxwell's equations. 5 
These equations are written as four coupled simultaneous first-order 
differential equations. The effective index N e for the single propagating 
HE n mode is found by solving the characteristic equation derived 
from the matching conditions of the field components at the core 
cladding interface. The dispersion is then calculated from N e . In 
addition, we find the conditions under which the TEoi and TAfoi modes 
propagate. 

Our computing procedures do not impose any restrictions on the 
index profile of the fiber. 5,6 However, for the sake of simplicity and 
familiarity we choose the a-index profile. To confirm the analysis, we 
compare the computed results with known experiments when possible. 
The calculated results include the dispersion, the optimum core radius 
(defined as the value at which the zero total dispersion will occur for 
a given relative index difference, wavelength, and a value), and the 
corresponding cutoff frequency. 

To keep the accuracy as high as possible, material dispersion is 
included in the calculation from the onset. The subtle interaction 
between material properties and waveguide properties determine the 
propagation characteristics of the lightguide. Details of the numerical 
procedure can be found in the appendix. 

II. DESIGN OF A DISPERSIONLESS SINGLE-MODE FIBER 

In the design of a single-mode fiber, the relative index difference A, 
the operating wavelength A, and the cladding index N2 are normally 
specified. The objective of the design is to make the total dispersion 
for the HEn mode equal to zero, the total dispersion D t being defined 
as the time spread of a narrow pulse per fiber length L per source 
spectral width dX. Qualitatively speaking, this can be accomplished by 
adjusting the radius and profile shape so that the waveguide dispersion 
exactly counterbalances the material dispersion. In addition, the design 
must be such that the TE, TM, and other higher-order HE and EH 
modes are not propagating. In the analysis shown in the appendix, the 
number of modes allowed to propagate in the fibert core are not 
restricted. However, the single-mode analysis is nothing but a partic- 
ular case of multimode propagation and it can be easily achieved 
simply by setting the angular mode number M = and 1 in the 
elements of matrix A [eq. (18) in the appendix] while fixing the radial 
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mode number Q = 1. When M = 0, the cutoff frequency for the 77? i 
and TMoi can be determined. Below this cutoff point, only a funda- 
mental mode (HEn) propagate in the fiber. All the modes other than 
M = are hybrid modes. Further details are explained below. 

When eqs. (22) and (23) are substituted into eq. (21), determinant 
(21) decomposes into two uncoupled equations, giving the TMand TE 
modes, 

TM mode: 

(j8?A2i - K(fi)YoAn) = . ^ 

(0?A 22 - KttihbAu) 



TE mode: 

(/8? A32 - yoA4 2 ) 



= 1. (2) 



(jBfAai - 70A41) 

When N\ is the maximum value of the index of refraction in the core 
and Nt is the cladding index, the effective index of the bounded modes 
must satisfy the condition N 2 <N e <Ni. Therefore, from eqs. (1) and 
(2) the effective index N e (0, 1) for the TMoi and TE i mode can be 
found. 5 The notation N e (M, Q) refers to the effective index of a mode 
having angular mode number M and radial mode number Q. 

In practice, the core radius is changed until N e (0, 1) is equal to N 2 . 
This gives us the cutoff radius a c - for the TE01 and TMoi modes. The 
normalized cutoff frequency is then calculated from 



1ma c 



Vc-Zp-JNl-Nl (3) 

A 

Now if V < V c , only a single mode (HEn) propagates. The effective 
index N e for this mode is calculated in the single-mode region as a 
function of the radius and is the principal quantity used to determine 
the dispersion of the HE] 1 mode. 

The relationship between the effective index N e (l, 1) and the group 
index N e is given by 

, dN e ... 

N« = N e -A— . (4) 

The propagation time of a pulse through a fiber of length L and group 
index N g is 

L 

C 

If there is a variation of N e with wavelength, there is a dispersion or 
spread in pulse width. Thus 
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t-=:N*. (5) 



where d\ is the spectral width of the source. 
We can substitute N g from eq. (4) into eq. (6) and get 

*— Jw^. (7) 



The total dispersion D t is 



_ dt 1 



Thus, 



* d'N m 

D ' = -cHxf (9) 



Obviously, we need to solve Maxwell's equations for N e and then 
calculate d 2 N e /d\ 2 to obtain Dt. 

Following the procedure described in our earlier papers, it is con- 
venient to describe the dispersive properties of the cladding by a 
modified Sellmeier formula. 5,7 Thus, 

where / = 0.035. The coefficients C, are given in Table I. Figure 1 
shows the variation of the silica cladding's refractive index, for the 
lightguides considered here. 

Since the relative index difference A is rather small, typically ranging 
from 0.2 percent to 0.8 percent, the core center index M can be written 
as 

N 2 

To obtain a feel for the size of the numbers involved, the following 
data are useful: 

A = 1.33 /an, 

N 2 = 1.446925, 

M = 1.449825 [as calculated by eq. (11)]. 



Table I 



Co = 1.4508554 
Ci - -0.0031268 
C 2 = -0.0000381 
C 3 = 0.0030270 
C 4 = -0.0000779 
C 6 = 0.0000018 
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Fig. 1 — Refractive index of silica as a function of wavelength. 



Since we assume a wavelength dependence for N2, described by eq. 
(10), Ni as calculated by eq. (11) is also wavelength dependent. The 
following profile formula is particularly useful and will be employed in 
the fiber design: 



N(r) = Ni[l - Ar a ]. 



(12) 



A few comments seem in order: When a — > 00 the profile is rectangular, 
when a = 2 it is parabolic, when a = 1 it becomes linear, and when 
a < 1 the profile develops a cusp. Figure 2 illustrates these shapes. 

In practice, one specifies A, X, and the profile parameter a. The 
computer then searches for a core radius that makes D t = by means 
of Mueller's iterative method. If it can be assumed that the total 
dispersion is the sum of waveguide dispersion D w (\, a, a, A) and 
material dispersion D m (X), then the physics of the problem is easy to 
understand. As will be shown later, this separation is quite accurate. 
Thus, a dispersionless single-mode fiber is one in which material 
dispersion is exactly balanced by waveguide dispersion. This is illus- 
trated by Fig. 3. 

Curve D m (\) shows the material dispersion as a function of wave- 
length for a single-mode lightguide with A = 0.2 percent and a rectan- 
gular profile. Note that it changes sign and passes through zero at 
about 1.27 /an. Also plotted are waveguide dispersion curves for three 
different radius values. Note that all three curves are negative and that 
as the radius decreases the magnitude of the dispersion increases at a 
given wavelength. 

From Fig. 3 we can draw a number of conclusions. First, the 
lightguide with the materials properties shown can be made totally 
dispersion free only at wavelengths longer than 1.27 fan. Second, as 
the wavelength gets longer the guide radius must get smaller. Finally, 
at longer wavelengths a much larger amount of material dispersion 
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Fig. 2 — Index profile shapes for various power laws. When a — * oo the profile is 
rectangular, when a — 2 it is parabolic, and when a — 1 it is linear. If a is less than 1, the 
profile develops a cusp. 



must be compensated for by waveguide dispersion. This requires 
greater precision in the waveguide parameters than when the guide is 
designed to operate at the zero of material dispersion. Hence, it is 
desirable to have a variety of materials available for lightguide design 
purposes. 

Finally, a typical total dispersion curve D t , illustrating lightguide 
operation at 1.3 /im, is plotted in Fig. 3. 

III. SINGLE-MODE LIGHTGUIDES WITH a PROFILES FOR 1.33-jum 
OPERATION 

Single-mode lightguides that are designed to operate at 1.33 jum with 
A ■ 0.2 percent are of much current interest. The cutoff frequencies 
V c (a) for the a- index profiles can be found from eqs. (1) and (2). Figure 
4 plots the results. There is little change in V c until the profile 
parameter is below 10, then it increases rapidly. At a = 1, V c and V op t 
are 4.383 and 2.724, respectively. The data of Fig. 4 agree well with 
that in the literature. 8,9 

To determine the radius a opt for zero total dispersion, a computer 
search is done for values less than a c , where a c is the cutoff radius [see 
eq. (3)]. Figure 5 shows the result for three a values, namely 100, 2, 
and 1. The radii turn out to be 4.142, 5.725, and 6.294, respectively. 
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Note that the total dispersion D t is most sensitive to the radius for 
a = 100, and least sensitive when a — 1. For design purposes, we plot 
a op t as a function of the profile parameter a in Fig. 6. This plot also 
includes the cutoff wavelength A c , where A c is defined by the following 
formula: 



Ac-(Vopt/Vc)-X. 



(13) 



Note that both V op t and V c are employed in the calculation. Therefore, 
to allow only one mode to propagate in the core, the operating 
wavelength must be longer than the cutoff wavelength A c . For a step 
index profile with A = 0.2 percent and X = 1.33 /mi, A c is close to 1 /xm. 
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Fig. 3 — A plot showing material dispersion D m , waveguide dispersion D u , and total 
dispersion D,. Note that because of the difference in sign between the material dispersion 
and the waveguide dispersion it is possible for one to cancel the other at a particular 
wavelength. 
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Fig. 4 — A plot of V e and V op t as a function of the profile parameter. 

The calculations described thus far in this section have made no 
assumption concerning the separation of the total dispersion D t into 
its component parts D m and D w . Our numerical procedure simply 
searches for a zero in total dispersion and does not assume that Dt = 
D m + D w . 

However, the material dispersion contribution can be calculated 
directly by means of eqs. (10) and (9), where N 2 is substituted for N e . 
The waveguide dispersion can be found by solving Maxwell's equations 
numerically via the procedure outlined in the appendix, with N2 being 
wavelength independent. 

Figures 7a through 7c show the results for a = 100, 2, and 1. In all 
three figures the material dispersion, which is of course independent 
of core radius, is represented by the horizontal dashed line. The 
negatives of the waveguide dispersions are dependent on core radius 
and are the solid lines in the figures. The intersection of the solid and 
dashed lines yields the optimum radii. As is evident from the figures, 
the radii are identical with the previous calculations. 

IV. A COMPARISON WITH EXPERIMENTAL DATA 

Miya, Terunuma, and Hosaka 10 have fabricated Ge(>2 doped single- 
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mode lightguides for the 1.3 /on and 1.5 fim region. Their experimental 
data are very useful to check the theory presented in this paper. They 
give data on a fiber with A = 0.2 percent and a zero total dispersion at 
X = 1.33 jum. A second fiber has A = 0.74 percent and a zero total 
dispersion at A = 1.54 /mi. Both are step-index fibers. Their experi- 
mental data are measured with a fiber Raman laser excited by a Q 
switched, mode-locked Nd:YAG laser. The spectral range studied is 
1.1 to 1.7 fim. The optimum radius is 4.142 /xm for A = 0.2 percent and 
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a. CORE RADIUS IN MICROMETERS 

Fig. 5— Total dispersion versus core radius for lightguides with a = 1, a = 2, and 
a* 100. Note that the total dispersion curve is steepest for a = 100 and shallowest for 
«-l. 
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Q. PROFILE EXPONENT 

Fig. 6 — A plot of cutoff wavelength \ c and optimum radius a opt versus a. This data is 
very useful for design purposes. 

A = 1.33 /mi. The optimum radius is 2.30 /mi for A = 0.75 percent and 
A = 1.53 /mi. Figure 8 compares their experimental results (circles and 
triangles) with our calculations (solid and dashed lines). The agreement 
is excellent. 

V. DISCUSSION 

At a wavelength of 1.33 /mi, dispersion-free lightguides can be 
designed with improved properties over the rectangular profile guides. 
Our calculated results of a op t and cutoff frequency V c (a) both show a 
rapid increase as the profile parameter a drops below 10. For a 
triangular index profile (a = 1), the core size is over 50 percent larger 
than the size of a step-index core. Thus, due to ease of handling, 
coupling to the source, and splicing, a larger core size is desirable for 
the same value of A and A. 

As the profile parameter a drops below 10, V op t also begins to 
increase substantially. It is significant that in the case of a = 1, V opt ~ 
2.7, which is a substantial increase over the step-index case, Vopt ~ 
1.8. 1U2 

In addition, Fig. 5 shows that a lightguide with a triangular profile 
is less susceptible to variations in core radius than a step-index guide. 
Therefore, because of the larger core size and less sensitivity in its 
tolerance, a triangular profile may be less difficult to fabricate than a 
rectangular profile. 
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APPENDIX 

To study the propagating modes in a fiber, we first assume that the 
cross section of the fiber is concentric in the core and cladding and 
that the core diameter 2a is much smaller than the outer diameter D 
(Fig. 9). Further, we assume that the core is a lossless, radially 
inhomogeneous medium with a scalar permittivity e. A relative per- 
mittivity (dielectric constant) is defined by k = €/c , where e is the 
value in free space. However, as is well known in a nonconducting 
medium, the permittivity becomes equal to the square of the index of 
refraction N. For source-free fields in a dielectric waveguide, Maxwell's 
equations reduce to 

V X E - -no — , 

(14) 



dt 



dE 

V X H = eofc — , 
at 
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Fig. 7(a) — A plot of waveguide dispersion and material dispersion as a function of 
core radius, a = 100. 
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Fig. 7(b) — A plot of waveguide dispersion and material dispersion as a function of 
core radius, a = 2. 



where E(R, t) and H(R, t) represent the electric and magnetic fields. 
fio is a scalar permeability in free space. 

In a cylindrical coordinate system R = {R, <$>, z) , the solution of eq. 
(14) is described by the vector components of two fields, {Er, E+, E 2 ) 
and {Hr, Hq, H z ). Among these components we are primarily inter- 
ested in finding the tangential components {E^, E 2 ) and {H+, H z ). 
These will be continuous through the core-cladding interface. Once 
the tangential components are known, the radial components Er and 
Hr can be obtained from these components. To find a complete set of 
bounded modes we consider the following form of solution: 13 



E(p, <J>, 2, t) = E(p) exp[iM - M<f>- /3z)], 



(15) 



where w and M are angular frequency and mode number, and ft is 
propagation constant along the z axis. 

Introducing a new variable A, the field components are written as 
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A = 



A, 




E z 


A 2 




-iZ pH+ 


A3 




pE* 


Aa 




-iZoH z 



(16) 



where p = KR — 2irR/\, Zo is the wave impedance defined by (/ioAo), 1 ' 2 
and A is a wavelength. 

Substituting eq. (16) into eq. (15) yields a set of the first-order 
ordinary differential equations. 1314 In the following equation, 



dp p 



(17) 



the operator A (p) is a 4 X 4 matrix that contains important property 
constants and parameters, for instance, the index of refraction N, 
angular mode number M, and the effective index N e defined by fi/k. 
The 4x4 matrix A (p) can be written as 



/ 
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Fig. 7(c) — A plot of waveguide dispersion and material dispersion as a function of 
core radius, a = 1. 
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Fig. 8 — A comparison of theory (dashed and solid lines) with experiment (see text). 



(18) 



The method of numerical computation has been given in our earlier 
work, 5 which briefly gives the procedure for finding the solutions and 
emphasizes what is needed to describe a single-mode fiber design. 

Equation (17) is a familiar type of differential equation. A fourth- 
order Runge-Kutta method is used to solve it numerically. A concise 
description is as follows. First, we know that there are two possible 
solutions in both the core and cladding, and they are linearly inde- 
pendent. Therefore, a general solution will be a linear combination of 
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two solutions, denoted by A, in the core and I\ in the cladding. 
In the core region, 

A, = AiKi + A 2 A i2 , (19) 

and in the cladding region, 

T, = AsAia + A 4 A,- 4 . (20) 

The solutions Aa and A l2 are transferred to the core-cladding inter- 
face through the operation of eq. (17), starting with two initial condi- 
tions given at the center of the core. ' 

At the interface p = p„ each has to be matched with the solution in 
the cladding to satisfy the continuity condition mentioned earlier. 
Thus, the matching condition yields a set of homogeneous equations 
for the constants At. Hence, the characteristic equation is the vanishing 
determinant of the system of equations, 



det|A,y|=0. 



(21) 




Fig. 9— Cross section of a fiber and cylindrical coordinate system. R or p is the radial 
coordinate, <f> is the angle, and z is the axial direction (p = KR ). 
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Particularly, considering M = 0, we can write the solutions in the core 
and cladding as follows. The two core solution vectors at p = p, are 
An and Aaipd- Expressed in terms of their components, they are 



Aa = 



A„ 
A 2 i 
A31 
A41 



and 



A l2 = 



A12 
A22 
A32 
A 42 



(22) 



And the two cladding solutions at p = p, are taken as 



Ai3 = 



01 

k(M) 






WoiSi) and A l4 = 







Yo(fi) 

A 



Wo(&), (23) 



where 



01 = [N 2 e - K($i)] 



1/2 



ft = PiPi, 



yo = ft-[Jr&(fc)/J£b(ft)], and W (ft) - &>(&)/$. 

The prime in the above equation denotes differentiation with respect 
tof. 
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